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Abstract: 

The Harmonic Einstein equation is the vacuum Einstein equation supplemented by a gauge 
fixing term which we take to be that of DeTurck. For static black holes analytically continued 
to Riemannian manifolds without boundary at the horizon this equation has previously been 
shown to be elliptic, and Ricci flow and Newton's method provide good numerical algorithms 
to solve it. Here we extend these techniques to the arbitrary cohomogeneity stationary case 
which must be treated in Lorentzian signature. For stationary spacetimes with globally 
timelike Killing vector the Harmonic Einstein equation is elliptic. In the presence of horizons 
and ergo-regions it is less obviously so. Motivated by the Rigidity theorem we study a class 
of stationary black hole spacetimes, considered previously by Harmark, general enough to 
include the asymptotically flat case in higher dimensions. We argue the Harmonic Einstein 
equation consistently truncates to this class of spacetimes giving an elliptic problem. The 
Killing horizons and axes of rotational symmetry are boundaries for this problem and we 
determine boundary conditions there. As a simple example we numerically construct 4D 
rotating black holes in a cavity using Anderson's boundary conditions. We demonstrate both 
Newton's method and Ricci flow to find these Lorentzian solutions. 
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1. Introduction 

Stationary vacuum black holes in spacetime dimension D > 4 have received much attention 
over the past decade, and many remarkable results have emerged [1, 2, 3]. In asymptotically 
flat D = 4 spacetime Kerr is unique. In D > 4 beyond the Myers- Perry solution [4] which 
generalizes Kerr, the space of solutions thought to exist is very exotic [5] and there may yet 
be further solutions discovered. String theory famously requires more than 4 dimensions for 
consistency, and hence provides a powerful physical motivation to think about these exotic 
solutions. Perhaps the strongest motivation to consider higher dimensions and also exotic 
boundary conditions is the AdS-CFT correspondence, which relates many of these interesting 
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solutions and other gravitational phenomena to physical phenomena found in certain strongly 
coupled conformal field theories and gauge theories where direct field theory calculations are 
extremely hard. The majority of the solutions conjectured to exist are unlikely to be found 
analytically using current techniques. The problem of solving the Einstein equations to find 
them reduces to solving partial differential equations for the components of the metric. Since 
these exotic solutions are of cohomogeneity greater than one it is not clear this can be done 
analytically in general. Thus we are forced to consider numerical methods to construct such 
solutions. The aim of this paper is to find numerical methods that are geometrically elegant, 
apply for arbitrary cohomogeneity and topology, and that function well in practice. 

From the classic AD stationary uniqueness theorems (see for example [6]) we understand 
that the stationary problem of finding vacuum black holes should be thought of as an elliptic 
PDE system which is treated as a boundary value problem, where asymptotic conditions 
are imposed together with regularity at the horizon. Rigidity implies the existence of ax- 
isymmetry, and then the problem reduces to one depending non-trivially on two coordinates 
('cohomogeneity two'). Taking the Weyl-Papapetrou form for the metric then manifests the 
PDE problem as an elliptic one. Black hole uniqueness theorems have been extended to D > 4 
in the stationary case of D — 2 commuting Killing directions [7, 8] where by reducing on these 
Killing directions the problem can be viewed as a 2D elliptic one on the orbit space of these 
isometries. For non-trivial dependence on only two dimensions the stationary numerical prob- 
lem has been extensively studied in the context of relativistic stars for decades [9, 10, 11, 12], 
using formulations based on the Weyl-Papapetrou form of the metric where a sufficient subset 
of the Einstein equations are elliptic. These methods have been applied to exotic static and 
stationary AD black holes [13, 14]. In [15, 16] it was shown how to apply cohomogeneity two 
methods to higher dimensions, without having {D — 2) commuting Killing vectors. A form for 
the metric analogous to the Weyl ansatz is taken and the Einstein equations reduce to a set 
of elliptic ones for the metric components, together with constraints which may be satisfied 
provided appropriate boundary conditions are imposed. 

The main challenge to formulating the problem for general cohomogeneity is to present 
the Einstein equations as an elliptic system in an elegant and covariant manner. In this 
paper we build on the previous algorithms proposed in [17] to find static vacuum black holes. 
In this method a static solution is analytically continued to a Riemannian manifold. If we 
are considering a black hole, then provided it has a single non-extremal horizon or multiple 
horizons sharing the same (non-zero) surface gravity, then the periodicity of Euclidean time 
can be chosen such that this Euclidean signature manifold has no boundary at the location of 
the horizon and is perfectly smooth there. The only boundaries in the problem arise from the 
region far from the horizon where one may consider various boundary or asymptotic conditions 
[18]. An elegant feature of this approach is that static vacuum Lorentzian geometries can then 
be treated simply as a particular case of the problem of finding Riemannian Einstein metrics 
which is more general and has other interesting applications [19, 20]. Instead of solving the 
Einstein equation directly, one solves the Harmonic Einstein equation with DeTurck's choice 
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of gauge fixing term. 1 For Euclidean signature this is elliptic and can be treated as a boundary 
value problem. The Ricci-DeTurck flow or Newton method provide explicit algorithms to find 
solutions [17]. Numerical Ricci flows have been considered in [21, 22, 23, 20, 24, 17]. 

This approach is closely related to the use of generalized harmonic coordinates in the 
hyperbolic context [25, 26]. One solves the Harmonic Einstein equation to yield a solution 
to the Einstein equation in generalized harmonic gauge. There is an important and subtle 
difference between the hyperbolic and elliptic contexts. For dynamics the vanishing of this 
gauge fixing term is imposed on the initial data, and the contracted Bianchi identity implies 
it remains zero at later times. Whilst one solves the Harmonic Einstein equation, one only 
considers this equation in the neighbourhood of a solution to the Einstein equation. In 
contrast in the Riemannian elliptic case one does not have this luxury. One must impose 
physical boundary conditions and then solve the Harmonic Einstein equation subject to these. 
Since one does not know the solutions which correspond to solutions of the Einstein equation, 
the whole problem being to find these, one cannot consider the Harmonic Einstein equation 
only in the neighbourhood of such Einstein solutions, but must consider it generally subject to 
the boundary data. In principle there may exist solutions to the Harmonic Einstein equation, 
termed Ricci solitons, where the solution is not Einstein, and the DeTurck gauge fixing term 
does not vanish. In practice one can check whether the solution found is a soliton or not. 
However, the existence of solitons is rather constrained. Bourguignon [27] has proven on a 
compact manifold without boundary such solitons cannot exist, and in [18] it was argued using 
a simple maximum principle that in this static case with various boundary and asymptotic 
conditions no solitons exist. 

Our aim here is to develop general methods to find stationary black holes. However this 
Riemannian approach for static solutions cannot be applied to the stationary case since there 
is no analytic continuation of the stationary vacuum Einstein equations. The key observa- 
tion is that the Harmonic Einstein equation of the Riemannian signature black hole gives 
precisely the same equations for the metric components if one continues back to Lorentzian 
signature provided one uses coordinates adapted to the static symmetry. Whilst in general 
the Lorentzian Harmonic Einstein equation has hyperbolic signature, this shows that when 
one considers the Harmonic Einstein equation restricted to static Lorenztian metrics it is 
elliptic and can be solved as a boundary value problem. The only difference is that now 
the Lorentzian horizon represents a physical boundary - in the Riemannian continuation the 
static adapted coordinate chart has a boundary at the horizon, in analogy to the origin of 
polar coordinates. One may determine the boundary conditions in these static coordinates by 
considering smoothness in a chart that does cover the horizon (analogous to Cartesian coordi- 
nates). The key result of this paper is to show that the Harmonic Einstein equation restricted 
to a class of Lorentzian stationary metrics motivated by the Rigidity theorem [28, 29], general 
enough to consider wide classes of stationary black holes which include the asymptotically 

1 Note that in the previous work [17, 18] the 'Harmonic Einstein' equation with DeTurck gauge fixing was 
referred to as the 'Einstein-DeTurck' equation. We will use the 'Harmonic Einstein' terminology here to avoid 
potential confusion with the 'Ricci-DeTurck' flow. 



- 3 - 



flat case, is elliptic. The horizon must be thought of as a boundary, and there are well posed 
boundary conditions that can be given there. The resulting elliptic system can be solved by 
Ricci flow on this space of Lorentzian stationary metrics, or by Newton's method. 

For simplicity we restrict our discussion here to vacuum solutions with no cosmological 
term and only non-extremal Killing horizons. The structure of the paper is then as follows. 
We begin in section §2 by reviewing the approach to finding static solutions using the Eu- 
clidean continuation and Harmonic Einstein equation. Then in section §3 we show how to 
reinterpret this from an entirely Lorentzian perspective, demonstrating the equation is elliptic 
and that regularity conditions can be given at the horizon, which in Lorentzian signature ap- 
pears as a boundary. In §4 we consider the character of the Lorentzian stationary Harmonic 
Einstein equation. We show that for globally timelike stationary Killing vector (i.e. station- 
ary spacetimes without horizon or ergo-regions) this is again elliptic. From this analysis it is 
clear that the threat to ellipticity in the stationary case is due to ergo-regions. In §5 we make 
an ansatz for stationary black hole spacetimes motivated by the Rigidity theorem, general 
enough to include the asymptotically flat case. This is the ansatz that Harmark [30] has used 
to consider classification of higher dimensional black holes. Following from Rigidity, Killing 
symmetry in the directions of angular or linear motion of the horizon is assumed, and this 
is crucial for maintaining ellipticity in the presence of ergo-regions. Based on the analytic 
treatment of the stationary problem for the uniqueness theorems, we assume the orbit space 
is a smooth Riemannian manifold, and the horizon and any axes of symmetries of the angular 
motion are boundaries of this space. We are then able to conclude that the Harmonic Einstein 
equation is elliptic when restricted to such stationary spacetimes. We determine the boundary 
conditions at the horizon and rotational symmetry axes in an analogous manner to those in 
the Lorentzian static case, and show they are compatible with finding non-soliton solutions 
to the Harmonic Einstein equations ie. solutions of the stationary Einstein equations. We 
conclude the paper in section 6 by giving a simple numerical example of these techniques. We 
numerically find AD rotating black holes in a cavity. We use Anderson's boundary condition, 
taking the boundary to be conformal to the product of time and a round sphere with fixed 
trace of its extrinsic curvature [31]. We demonstrate the use of both Newton's method and 
Ricci flow as algorithms to find the solutions. 

2. Review of the Harmonic Einstein equation, Ricci flow and Newton's 
method for static geometries 

We now review in more detail the framework for finding static vacuum black holes that 
we wish to extend to the stationary case. Since we will consider vacuum solutions with 
no cosmological term we are interested in finding Ricci flat static metrics. We analytically 
continue time to give a Euclidean metric. We take all the Killing horizons to have the same 
surface gravity so that we may choose Euclidean time to have the appropriate period such 
that we have a smooth Riemannian manifold, and no boundaries associated to the horizons. 
Thus the problem is now to find Ricci flat Riemannian metrics with a U(l) isometry which 
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is generated by a hypersurface orthogonal Killing vector, and which has fixed points of the 
isometry at the locations which continue to Killing horizons. 

Following our previous work, instead of solving the Ricci flatness condition R^ u = 0, we 
instead introduce the tensor, 

which we will refer to as the Ricci-DeTurck tensor, where the vector field £^ is constructed 

as, 

e = ff* (TJ - Tj) (2.2) 

from a fixed smooth reference connection rjj 1 on the manifold. The term V(^£j,) is the 
'DeTurck term' introduced by DeTurck in the context of Riemannian geometry and later 
used by him to show Ricci flow is parabolic [32]. We note that the vector field £ is globally 
defined, being the difference of our Levi-Civita and reference connection. This term is a global 
analog of the local term considered in generalized harmonic gauge fixing [25, 26] which is now 
extensively used in numerical simulation following Garfinkle's work. 

We focus on solving the 'Harmonic Einstein equation' R^ u = 0. The Ricci-DeTurck 
tensor has the property that it exhibits a simple structure in its second derivative terms, so 
that, 

K> = R »» ~ V (f&) ~ -\a aP d a d p g^ + ... (2.3) 

where . . . refer to terms of lower order in derivatives. Then for a Riemannian signature 
manifold the Harmonic Einstein equation is clearly elliptic, and can be solved as an elliptic 
boundary value problem in the metric components g^ v . Whilst the Ricci tensor R^ u shares 
the static symmetry of the metric g^ u , R^ u does not unless we choose the reference connection 
appropriately. We do so by taking the reference connection to be the Levi-Civita connection 
of a smooth reference metric, g^ u , so that, 

= sT (vtfua ~ ^V Q (7^J (2.4) 

and we take this reference metric to share the static isometry. Here, V is the covariant 
derivative with respect to the metric g^y. This ensures R^ u is symmetric with respect to the 
static isometry. Both metric and reference metric continue to Lorentzian spacetimes with the 
same static Killing horizons with equal surface gravities with respect to Killing time. 

The aim is to solve the Harmonic Einstein equation to yield a Ricci flat solution with 
^ = which specifies that the metric is presented in the generalized harmonic gauge defined 
by our reference metric. However in general there will exist solutions which are not Ricci flat, 
but solve R^y = V^j^j,) so that £ M does not vanish. Such solutions are called Ricci solitons. 
Fortunately the existence of Ricci solitons is quite constrained as we shall now discuss. The 
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contracted Bianchi identity for the Harmonic Einstein equation implies £ obeys the elliptic 
linear equation, 



(2.5) 



The boundary conditions for the metric may be taken as giving boundary conditions for the 
vector field £. Let us consider a vector field with the same boundary behaviour as that of 
£ At . We expect that our boundary conditions are such that the problem, 



is a well posed elliptic linear PDE in the vector components v^, and that = is a solution. 
For example, taking asymptotically flat boundary conditions then £^ — > in the asymptotic 
region [18]. Alternatively taking Anderson's boundary conditions where the conformal class 
of the metric is fixed on a boundary, together with the trace of the extrinsic curvature [31, 18] 
then one imposes £^ = 0. For the boundary conditions relevant for a vacuum Randall- 
Sundrum brane [33, 34], the normal component of £ M to the brane boundary vanishes, and 
the tangential components of £^ have Neumann boundary conditions. In all these three 
examples the associated problem has elliptic boundary conditions for the vector v^. In the 
first two these ensure vanishes, whereas in the latter one, they are consistent with 
vanishing but do not impose it. 

A necessary condition for a Ricci soliton to exist is that the linear elliptic operator O 
must have a non-trivial kernel. For given boundary conditions on the metric we may view this 
as an obstruction to the existence of Ricci solitons. A well known result is that for compact 
manifolds without boundary there are no Ricci soliton solutions for any choice of vector field. 
Using a simple maximum principle based on the contracted Bianchi identity one can show 
that with a variety of boundary conditions, including asymptotically flat, Kaluza-Klein and 
Anderson's boundary conditions, then again no Ricci solitons can exist [18]. 

However, even assuming Ricci solitons do exist, then the well-posedness of the Harmonic 
Einstein equation implies that generically 2 solutions should be locally unique. Hence a Ricci 
flat solution cannot be perturbatively close to a soliton solution. Therefore even if solitons 
exist, it is straightforward in principle to identify whether a solution that has been found is 
a soliton or not, the most obvious way being to compute £ and determine if it vanishes. 

Here we will now briefly review Anderson's boundary conditions since we will use them 
later in our numerical example. Anderson has shown the surprising result that one cannot 
impose the most obvious Dirichlet boundary condition, namely fixing the induced metric on 
a boundary, since it does not give rise to a regular elliptic system. Instead he has argued that 
one can fix the conformal class of the induced metric on a boundary together with the trace 

2 There may be special points in moduli space of Ricci flat solutions where a normalizable zero mode of the 
linearisation of the Harmonic Einstein equation exists, and this lifts to the non-linear equation to generate a 
deformation to a soliton solution. If a branch of solitons meets a branch of Ricci fiat solutions this would be 
the case. 




= 



(2.6) 
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of the extrinsic curvature [31]. The metric in the vicinity of such a boundary can be written 

as, 

ds 2 = Ndr 2 + N a drdx a + g ab dx a dx b (2.7) 

where the boundary is the hypersurface r = and the remaining tangential coordinates are x a 
for a = 1, . . . ,D — 1. For a static non-extremal black hole continued to Euclidean signature, 
the boundary will be a product of the Euclidean time circle with some spatial geometry. 
Fixing the conformal class of the induced metric and trace of the extrinsic curvature locally 
gives D(D — l)/2 conditions on the metric functions. However since the metric has D(D + \) /2 
independent components we require another D conditions. Recall our previous discussion, 
that on any boundaries or asymptotic regions we must ensure that £ has a behaviour such 
that the associated linear problem in equation (2.6) defined above is well posed and has trivial 
solution = 0. In order to achieve this the remaining D conditions are fixed by imposing 
^ = on the boundary. 

For Riemannian metrics we may solve the Harmonic Einstein equation using two distinct 
methods. Firstly we may use the DeTurck flow, 

^ v = -2R% (2.8) 

beginning with an initial guess metric, and then flowing in the auxiliary time A, hoping to 
end at a fixed point of the flow, = 0. The flow is strongly parabolic provided R^ is 
elliptic. Furthermore for our choice of reference metric it preserves the static isometry of 
g. A very elegant feature of this flow is that it is diffeomorphic to the Ricci flow ^g^v = 
— 2R lll/ , since the DeTurck term V^^) just introduces an infinitessimal diffeomorphism at 
each point along the flow. Hence while the flow in the space of metrics will depend explicitly 
on the choice of reference connection T, geometrically the flows do not depend on this choice. 
An important consequence is that the basin of attraction of a fixed point of the flow is 
geometrically invariant, and doesn't depend on the choice of T. One important drawback of 
this method is that, as discussed in [22] vacuum black holes of interest are often unstable fixed 
points of Ricci flow, due to the existence of Euclidean negative modes of the Lichnerowicz 
operator [35]. However, the Ricci-DeTurck flow may still be used as an algorithm. If the fixed 
point of interest has n negative modes, then as described in [22], an n parameter family of 
initial guesses must be taken, and the n parameters tuned in the sense of a shooting problem 
in order to reach the fixed point of interest. 

The second method is to consider solving the equations Rr[ v in a non-local way, using the 
canonical method to solve non-linear equations, the Newton's method. Writing the linearisa- 
tion of R H as, 

R%[g + e5g](x) = R%[g](x) + eA[ 9 ](x)^ 9a/3 (x) + 0(e 2 ) (2.9) 

we may think of this as an (infinite dimensional) vector equation, with O(g) a matrix, with 
indices constructed from the spacetime index pair [iv and the point x on the manifold. In 
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practice one will introduce a truncation of the spacetime, for example using finite difference 
or spectral methods. Let us denote the finite dimensional collective index which represents 
the spacetime and position/mode using A, B, so that we may write g^ u (x) as the vector g^- 
Then we write the above equation as, 

Ra (9 + eSg) = R% (g) + eA(g)£5g B + 0(e 2 ) (2.10) 

(i) 

The Newton method iteratively improves a guess metric g A as, 

^ 1) -»2 ) = -(A( & W)- l )/«fa/ (i) ) 

As for the Ricci flow method, the Newton method preserves the static isometry of the metric. 
This method has the important advantage over the Ricci-DeTurck flow method that it is not 
sensitive to negative modes of the Lichnerowicz operator. However it does assume that the 
linear problem A ■ v = R H can be solved for v. In practice robust methods exist to solve 
such (finite dimensional) linear systems, such as biconjugate gradient, which are insensitive 
to the spectrum of A, provided there are no zero modes. Thus a single initial guess will 
suffice, rather than having to tune a family of initial guesses. An important disadvantage 
of the Newton method over the Ricci-DeTurck flow is that it is not geometric in the sense 
that the path taken by the algorithm in the space of geometries will depend explicitly on the 
choice of reference connection f . This implies that the basin of attraction of a solution will 
also depend on this choice. 

3. Static spacetimes from a Lorentzian perspective 

Instead of immediately now considering stationary spacetimes, it is instructive to first con- 
sider static spacetimes from a Lorentzian perspective. The Harmonic Einstein equation is 
not elliptic for a general Lorentzian manifold, and without ellipticity one would not expect 
to be able to impose the various boundary conditions that we require physically in a well 
posed manner. However, consider a chart away from any horizon which manifests the static 
symmetry, 

ds 2 = g^dx^dx" = -N(x)dt 2 + h ij (x)dx i dx j (3.1) 

so that N > 0. We may regard Euclidean time as being fibered over a base, the Riemannian 
manifold which we shall denote M with Euclidean metric ds 2 M = hij(x)dx % dxK Furthermore 
with the choice that our reference metric is also static with respect to d/dt, so that, 

ds 2 = g^dafdx" = -N{x)dt 2 + h ij (x)dx i dx j (3.2) 

again with N > and hij a smooth Euclidean metric, then Rjf v shares the static symmetry. 
Due to this static symmetry the Harmonic Einstein equation R^ v = thought of as PDEs for 
the metric components of g is invariant under an analytic continuation t — > r = it. Hence we 
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immediately see that the Harmonic Einstein equation restricted to Lorentzian static metrics 
and reference metrics is elliptic. Furthermore, Ricci-DeTurck flow yields precisely the same 
flow equations for the metric functions N and hij above in either Lorentzian signature, or 
under continued Euclidean signature. 

Explicitly the static Ricci-DeTurck tensor has components, 

+ ^N){d, t N) - \i k d k N - -Lh^^NMN) 



^ ~ ~ W2 (9iN)(9jN) - Ih^fi (^h^i)d m Nj (3.3) 

where indices are contracted and covariant derivatives V are with respect to the base metric 
hij. To avoid confusion we use the notation h 1 /^ for the inverse metric to hij, so that 

hikh^ 3 _^ = h^^hjk = 5 l k . The vector field £ l is the DeTurck vector of the base metric, 
namely, 

C = W k (f* ifc - f^) (3.4) 

where P - fe is the connection for hij and P - fc is the connection for hij. We see the terms that 
contain two derivatives acting on the metric components of g are of elliptic form provided 
h^ is a Riemannian metric. We also explicitly see that Rjf v is symmetric under the static 
isometry. 

In the exterior of any horizons we have N, N > 0. However now consider a non-extremal 
Killing horizon, where d/dt has fixed action. Now N and N vanish at the horizon. The base 
hij must remain a smooth Riemannian geometry where N = for the spatial horizon to have 
a well defined geometry, and we choose the same for hij. In the Riemannian case we know 
that if we have have a single non-extremal horizon, or multiple horizons with the same surface 
gravity, then we may make Euclidean time periodic as, r ~ r + 2ir/n for some appropriate 
choice of constant k, such that there is no boundary at the horizon(s) and the geometry is 
smooth there. Furthermore we know that Rjf v is a smooth tensor on this geometry, and hence 
the Ricci-DeTurck flow and Newton methods preserve the smoothness and lack of boundary 
at the Riemannian horizon. However, in the chart above with t — >■ r = it, the metric becomes, 

ds 2 = +N{x)dT 2 + h ij (x)dx i dx j (3.5) 

and the chart does not cover the horizon where N = 0. Such coordinates adapted to the static 
symmetry are analogous to polar coordinates, and fail at the polar origin, the horizon. To 
manifest the smoothness of the Riemannian manifold one must go to 'Cartesian' coordinates. 
Taking coordinates in the base adapted to the horizon such that x % = (r, x a ) where r = is 
the horizon, we write, 

ds 2 = +r 2 Vdr 2 + U{dr + r U a dx a ) 2 + h ab dx a dx b (3.6) 



R 



H 



rt t j 
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where the metric functions are functions of r and x a . Changing to coordinates, 



a = r sin kt 



b = r cos kt 



(3.7) 



provides a good chart covering the horizon, such that the metric components are smooth 
functions, provided that V, U, U a , h a b are smooth (C°°) functions of r 2 and x a , and, 



at the horizon r = [17, 18]. Precisely the same conditions will apply to the reference metric 
which is also required to be smooth with no boundary at r = 0. We now see that instead 
of using a good chart which does not manifest the static isometry, we might just as well use 
the original 'polar' chart (3.6) and simply treat the horizon as a boundary, and determine 
the boundary behaviour using the regular chart, namely that V,U,U a ,h a b are smooth in r 2 
and x a , and that V = k 2 U at r = where the constant k determines the angular period of 
Euclidean time, and hence the temperature of the solution. In practice, adapting coordinates 
to the static symmetry, and indeed any other isometries, is important numerically to simplify 
the problem, and this is precisely the approach taken in previous work [22, 17]. As noted 
above taking smooth coordinates we know preserves smoothness and lack of boundary 
at the horizon. Whilst in adapted coordinates we have a boundary at the horizon the same 
must be true, namely that the smoothness and regularity conditions above apply equally well 
to the tensor Rjf v which must be a smooth tensor on g and we may explicitly check that, 



where the functions f,g,g a an d r ab are smooth in r 2 ,x a , and in addition / = n 2 g at r = 0. 
However, the simpler way to see that this must be the case is to remember that in the smooth 
Cartesian coordinates (a, b, x a ) then g^ is smooth and hence Rj* v will be too, and since R^ v 
with our reference metric preserves the static isometry, then it follows that R^ must have 
the behaviour stated above. 

In Lorentzian signature, if we are to study only the exterior of the horizon, then the 
horizon should be regarded as a physical boundary. There is no analog chart to the Rie- 
mannian case where the boundary can be smoothly removed without introducing the black 
hole interior. However, since the Harmonic Einstein equation and the solutions for the metric 
components are invariant under continuation t —¥ t = it, then precisely the same boundary 
conditions for regularity apply in the Lorentzian case (3.1) as in the Euclidean case in the 
static adapted chart (3.5). Hence we may work directly in Lorentzian signature, where the 
equations are the same and so are elliptic, and then provide the same boundary conditions 
there for the metric components at the horizon, taking, 



V = k 2 U 




R H = + r 2 fdT 2 +g{dr + r g a dx a f + r ab dx a dx b 



(3.9) 



ds 2 = -r 2 Vdt 2 + U{dr + r U a dx a f + h ab dx a dx b 



(3.10) 



where r = at the horizon, V > outside the horizon and vanishes on it. Then at r = 
we again require V, U, U a , h ab are smooth in r 2 and x a , and that V = k 2 U. Now the constant 



k precisely gives the surface gravity of the Killing horizon with respect to d/dt. We may 
manifest the regularity of this horizon in a similar manner to the Euclidean case by performing 
a change of coordinates, 

a = rsinhftt, b = rcosh Kt (3-H) 

giving a chart with coordinates a, b, x a that now covers the t = slice of the Killing hori- 
zon, and whose metric components are smooth functions. The essential difference with the 
Euclidean case is that if we are interested in the exterior of the horizon, then the horizon 
remains a boundary in this good chart. 

Due to the invariance of the components of under t — >• r = it, the tensor shares 
the same regularity properties as the metric in the Lorentzian context too so that, 

R H = -r 2 fdt 2 + g{dr + r g a dx a f + r ab dx a dx b (3.12) 

near the chart boundary at r = where again f,g,g a and r ab are smooth in r 2 ,x a , and in 
addition / = K 2 g. In the Riemannian picture Ricci flow and the Newton method preserve 
smoothness and lack of boundary at the horizon. We see that equivalently in the Lorentzian 
picture we have the very nice property that Ricci flow and the Newton method will preserve 
the surface gravity of the solution. 

We conclude this discussion with some comments. One attractive feature of the Rie- 
mannian approach to static black holes where one removes the horizon boundary by taking 
periodic time is that there is no boundary associated to the horizon, and the boundary con- 
ditions imposed far from the horizon may be taken to fix the size of the time circle. Since the 
size of the time circle is interpreted as inverse temperature, we see that viewing the system 
as a boundary value problem we are naturally led to fix physical data. A very nice conse- 
quence of our static Lorentzian discussion is that although we must now view the horizon as 
a boundary, we are again naturally lead to impose physical data there, namely the surface 
gravity with respect to d/dt. Asymptotically or on a boundary away from the horizon we 
impose conditions to fix the value of the function N in (3.1), and this then determines the 
normalization of d/dt. Thus together these fix the physical data specifying the black hole, 
and moreover this data is preserved by the Ricci flow and Newton method. 

Since in practice (for example [17]) we choose to use adapted coordinates that manifest 
the static symmetry even when thinking about the Euclidean formulation of the problem, 
and then apply boundary conditions where the coordinates degenerate, one might think it 
makes no difference which signature we think about the static problem in. The mechanics of 
solving it will be identical in both. Whilst true, there is one important advantage to thinking 
about the static problem from the Lorentzian perspective, namely that one can consider 
multiple Killing horizons with respect to d/dt with different surface gravities, each of which 
is individually preserved by Ricci flow or Newton's method. Formally we could not previously 
consider this in the Riemannian CELS6, clS only one horizon boundary can be removed and made 
smooth by a choice of periodic time, the remaining horizons becoming conical singularities. 
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By employing the boundary conditions above in the static adapted coordinates of course we 
may now treat these conical singularities in the Euclidean context, but the motivation to 
consider Euclidean signature is somewhat diminished. 

4. Stationary vacuum spacetimes with globally timelike Killing vector 

We now wish to consider using the methods above to find stationary vacuum solutions. How- 
ever, now the spacetimes must be considered as Lorentzian from the outset as there is generally 
no real Euclidean section to a solution. In this section we begin by considering the case of 
stationary spacetimes with globally timelike Killing vector (i.e. no horizons or ergo-regions 
may exist) and will argue that the Harmonic Einstein equation is elliptic. Of course we are 
ultimately interested in black hole spacetimes which violate such a condition and in the fol- 
lowing section §5 we consider more general stationary spacetimes which allow horizons and 
ergo-regions. 

Consider the most general stationary metric with Killing vector T = d/dt, which we may 
write using coordinates adapted to the stationary isometry as, 

g = -N(x) (dt + A i {x)dx i ) + h ij (x)dx i dx j (4.1) 

Now under our assumption that T is globally timelike we have N > and we further assume 
that the function N is bounded. Physically this implies our spacetime has no Killing horizons, 
and also no ergo-region. Since det g^ v = —Ndethij we see that provided the metric g is 
Lorentzian and smooth, so that detc^ < and bounded, this implies that det hij > 0. We 
may then regard this metric as a smooth fibration of time over a base manifold A4 so that 
(A4 , h) is a smooth Riemannian manifold with Euclidean signature metric hij . 

Whilst one might imagine that for a stationary spacetime the natural way to think about 
the metric is using the ADM ansatz, it is worth noting that the metric above is not of this 
form, but rather should be thought of as a Kaluza-Klein ansatz over time. 3 Hence the base 
manifold (A4,h) is not the submanifold obtained by taking a constant time t slice of the 
spacetime. Rather it is the geometry one obtains by performing a Kaluza-Klein reduction in 
the time direction. As a result of this the second order derivative terms acting on the metric 
components g^ u in the stationary Harmonic Einstein equation then go as, 

R% ~ " \s^d a d pgiw + ... = —h^didjg^ + ... (4.2) 

and we see whilst the metric g^ u is indeed Lorenztian, since there is no dependence on the 
coordinate t, it is actually the metric hij that controls the character. This immediately implies 
that the Harmonic Einstein equation Rjf v = is elliptic since h is smooth and Euclidean 
signature. Thus the stationary problem reduces to an elliptic problem on the Riemannian 
base manifold M. 

3 Interestingly in the seemingly unrelated context of effective field theory used in the context of GR the 
utility of such a Kaluza-Klein ansatz has been emphasized over the ADM one [36]. 
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However, we also require that R^ v is a tensor that is symmetric with respect to the 
stationary isometry T. Without this, Ricci-DeTurck flow and Newton's method will not 
consistently truncate to the class of stationary metrics (4.1). In order that R~f v preserves the 
symmetry T, we choose the reference metric g to also be a smooth Lorentzian metric which 
is stationary with respect to the vector field T, so that, 

g = -N(x) (dt + A i (x)dx i ) 2 + h ij (x)dx i dx j (4.3) 

where we also assume here that T is globally timelike and bounded with respect to g so that 
N > and bounded. Then hij gives a second Riemannian metric on the same manifold 
Ai . Since R^ v preserves the stationary symmetry, the Ricci-DeTurck flow can be consistently 
truncated to a parabolic flow on the space of Lorentzian stationary metrics. Since this flow 
remains diffeomorphic to Ricci flow (subject at least to the normal component of £ vanishing 
on any boundaries), we arrive at the interesting result that we may apply Ricci flow to 
Lorentzian stationary spacetimes. Likewise the Newton method will preserve the symmetry. 

In a situation where the solution we wish to find has a stationary Killing vector that is 
globally timelike and bounded then nearby to that solution the character of the Harmonic 
Einstein equation will be elliptic. Subject to imposing suitable boundary conditions on any 
boundaries or asymptotic regions, one may use the Lorentzian stationary Ricci-DeTurck flow 
or Newton method to solve for the solution. One must start with an initial guess that 
has globally timelike bounded stationary Killing field T, and then provided that guess is 
sufficiently good the subsequent Ricci-DeTurck flow or Newton iterations preserve that T is 
globally timelike. 

We recall that for a solution to the Harmonic Einstein equation to be Ricci flat we require 
£ = 0. Thus we must ensure that our boundary or asymptotic conditions are compatible with 
vanishing £. Let us briefly consider the case of a Dirichlet boundary. The situation is entirely 
analogous to the static Euclidean case discussed above. Consider taking coordinates adapted 
to the boundary, so that, 

ds 2 = -N(x) {dt + A r (x)dr + A a (x)dx a ) 2 + Vdr 2 + V a drdx a + h ab (x)dx a dx b (4.4) 

so that x % = (r, x a ) and the boundary is at r = 0. Then fixing the induced metric would 
specify Dirichlet conditions for N, A a and h ab . Requiring that £*, £ r and £ a vanish then 
provides conditions for V, V a , and A r . Thus as in the static case we have precisely fixed 
both the induced metric and ^ = on the boundary. Following an analogous argument as 
in [18] we expect these boundary conditions are well posed for the elliptic Harmonic Einstein 
equation. Provided on other boundaries conditions are set consistent with £ vanishing we 
expect solutions where £ vanishes globally may be found (presuming they exist). Without 
the maximum principle of the static case [18] we cannot rule out solitons with non-vanishing 
£, so in a practical context one simply has to test a solution found to see if it is Ricci flat or 
a soliton. 

4 In the static case we have the inequality X7 2 (f> + £ M <9 M = V M £„V£" > where cj> — £| 2 . However in 
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For completeness we now explicitly give the stationary Ricci-DeTurck flow equations 
(which of course give the Harmonic Einstein equations at a fixed point d/dX = 0); 

=h km d k d m N+... 

+ ^~ h (-i)( d ^)(d k N) + h^qNiA* - A^FjmdkN A^h^^N^N) 



^ = V k F lk + Vi(A k e) + Vi(VM p - V p A p ) —FyPN + t*F u + — h^F^N 
- h k ^ 1} N(Ai - Ai)F ik F jm + h^Ap - A p ) 2 F ki d m N + V, ((^(^(An " An) 
+i(A" - i p ) + ^f! 1} (A m - A») 2 (Afc - A,)) d p iv) + Vi^iV^ - - A k )F 



jmj 



=h km d k d m h l3 + 

dhij ' „^ P 1 



^ = -21Jy + 2V (i £ i} +^(d i N)(d j N) - NFjF ki 

-hi k Vj( — ii,_ ] ,(./,,, ! -r //,/, \ , [ii ( _ | , 
1 



h ik V j{^ hk (-i)dmN) + h ik V j (M™ ) NF qrn (A<i - A*)) 



+^h ik V j (h^ 1) (AP - An 2 d m N) |(m j) 

where indices are contracted and covariant derivatives V are with respect to the base metric 
h{j, and we have defined the antisymmetric tensors F^j = diAj — djAi and Fij = diAj — <9jA- 
Again £ l is the DeTurck vector of the base metric defined as before in (3.4). 

We reiterate that the assumption the reference metric is also stationary with respect 
to T is responsible for R^ u preserving the stationary symmetry and hence the consistent 
truncation of the Ricci-DeTurck flow to the class (4.1) that we see above. However the 
further assumption that g has globally timelike bounded Killing vector T, so that N > and 
bounded, ensures that the equations are regular and there are no singular terms arising from 
vanishing or diverging N in the above. We note that the components of the DeTurck tensor 
R^ u may be derived from these DeTurck flow equations using the relations, 

9N - 2J? H 



dh 



ox 



~2{Rij + RuAiAj — RuAj — RjtAi) 



the stationary case this inequality does not appear to hold as in Lorentzian signature with only stationary 
symmetry V M £,/V M £" is not of definite sign. 
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5. Stationary black holes and the Harmonic Einstein equation 



We now proceed to consider the case of Ricci flat non-extremal stationary black holes. In the 
context of the discussion above now the norm of T will vanish either at the horizon itself, 
assuming that T is a globally timelike Killing vector (such as for certain Kerr-AdS black holes 
[37]), or outside the horizon at the boundary of the ergo- region. Since we are interested in 
the exterior of the horizon, in the first case we may treat the system described above for 
globally timelike T and now regard the horizon as a boundary of the base manifold Ai where 
suitable boundary conditions are required. However in the latter, more general case, outside 
the horizon but inside the ergo-region we have the norm of T > and hence det hij < 0. 
Now the base manifold in the previous section fails to be Riemannian and then our argument 
above that the Harmonic Einstein equation is elliptic fails. 

In order to make progress we must use the Rigidity property of stationary black holes, 
proved in D > 4 by Ishibashi, Hollands and Wald [28] and by Moncrief and Isenberg [29] for 
various asymptotics including asymptotic flatness. Assume there exists a stationary Killing 
vector T. Then the Rigidity theorem states that for a rotating non-extremal Killing horizon 
with topology ixE, for £ compact, there exists a Killing vector K that commutes with T, 
which is normal to the horizon. Furthermore there exist some number N > 1 of commuting 
Killing vectors R a , which also commute with T and generate closed orbits with period 2tt and 
K may be written in terms of these as, 



for some constants £l a . This result is physically significant as it ensures that the rotation of 
the horizon is generated by an isometry of the spacetime. Were this not the case one would 
expect gravitational radiation to be emitted from the region near the horizon and this would 
presumably violate the assumption of stationarity. 

Motivated by the Rigidity theorem we assume that our stationary spacetime, with sta- 
tionary Killing vector T, has iV additional Killing vectors R a for a = 1, . . . , N, which commute 
amongst themselves and with T and generate rotational or translational isometries. In the 
former case they generate closed orbits which we take to have period 2ir, and may have axes 
of rotational symmetry where R a vanishes. In the latter case they generate non-compact 
orbits. Let us take there to be a number of disconnected horizon components, H.\, . . . ,Hk- 
Then Rigidity implies that each component is a Killing horizon with Killing vector given by a 
linear combination of the isometries T and R a so that, Ku m = T+H,^^R a for some constants 
( w hich may be different for each component). 

As a consequence of these assumptions we may write the metric adapting coordinates to 
the isometries, so that using coordinates y A = {t,y a }, 



ds 2 = g^dXW = Gab(x) {dy A + Af(x)dx i ) (dy B + Af(x)dx j ) + h^dx^x^ (5.2) 



where T = d/dt and R a = d/dy a . In analogy with the stationary case in the previous section 
we see that the geometry may be thought of as a fibration of the Killing vector directions 



K = T + n a R, 



a 
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over a base manifold Ai with metric hij . This base manifold A4 is the orbit space of the full 
Lorentzian spacetime with respect to the isometries T,R a . If R a generates a compact orbit, 
then since we have chosen to normalise the period to 2ir, then the coordinate y a is periodic 
with y a ~ y a + 2n. Again we note that the metric hij is not that induced on a constant y A 
submanifold of the full spacetime, but rather is the metric one would obtain by performing a 
'Kaluza-Klein' reduction over these Killing directions. 

The full spacetime is Lorentzian and so exterior to the horizons det g^ u = det Gab det hij < 
0. On physical boundaries or asymptotic regions T is timelike, R a are spacelike and hence 
the fiber metric Gab is Lorentzian there, and consequently the base metric is Euclidean. At 
a horizon Ti m , since the norm of Ku m = T J r n^ im R a vanishes, then det Gab = 0. In addition 
at axes of symmetry associated to the fixed action of a compact R a , det Gab will vanish. 

Following the uniqueness theorem treatment of AD stationary black holes as an elliptic 
problem on the two dimensional Riemannian orbit space bounded by the horizon and axes of 
symmetry (as for example discussed in [6]) and its generalization to D dimensional metrics 
with (D — 2) commuting Killing vectors [7, 8] which is treated in the same manner, we make 
the follow key assumption; 

• We assume that the orbit space base manifold (M, h) is a smooth Riemannian manifold 
with boundaries given by the horizons and axes of symmetry of the R a that generate 
rotational isometries. 

A consequence of det hij > everywhere on M. (including the horizon and axis boundaries) , 
is that det Gab > everywhere on A4 with it vanishing only at the horizon or axis boundaries 
of A4. We note that detg MJ , = at the horizons and axes, as one would expect since the 
chart (5.2) breaks down there. As Harmark has discussed [30], the structure of M together 
with the data Qi, at the horizon boundaries, and the data of which R a vanishes at the axis 



boundaries defines a 'rod structure' for stationary spacetime and has been conjectured to 
classify higher dimensional black holes. 

5.1 Ellipticity 

We note that we have not considered the stationary Killing field T to be timelike. In the 
presence of horizons it will become null on the horizon or be spacelike if the horizon is 
surrounded by an ergo-region. We reiterate that in the previous section §4 it was precisely 
where T failed to be timelike that ellipticity would break down, since the base metric would 
fail to be Riemannian. The crucial observation is that for our class of stationary spacetimes 



■m 



(5.2) 



Rab = -\g aP d a dpg A B + ■■■ 



-\h mn d m d n G AB + ... 
l -h mn d m dnhij + ... 



(5.3) 
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where again the . . . represent lower than second order derivative terms. We see the equations 
have character determined solely by the metric hij, and by our assumption above that the 
base A4 is Riemannian, this is indeed elliptic. 

In analogy with the previous section 4, in order to ensure that R B U shares the symmetries 
of g we choose the reference metric g so that T, R a are again Killing with respect to it, and 
obey precisely the same assumptions as above for g. Thus we may write, 



ds 2 = g^dX^dX" = Gab(x) {dy A + Af(x)dx i ) (dy B + Af(x)dx j ) + h^dx^x^ (5.4) 



and we further assume that (M,h) is a smooth Riemannian manifold. Then the Ricci- 
DeTurck flow and Newton's method consistently truncate to the Lorentzian stationary space- 
times of the form (5.2). 

We must impose suitable boundary conditions. In addition to the boundaries that define 
the asymptotics we must also now treat the boundaries at the horizons and axes of symmetries 
of the rotational Killing vectors. We will shortly discuss these boundary conditions explicitly. 
Using the Ricci-DeTurck flow or Newton method if we start from initial data in our stationary 
class, then for small flow times or updates we expect to remain in this class. In particular 
we expect (A4,h) to remain a Riemannian manifold. Provided this condition holds for the 
solution of interest, and our initial guess is sufficiently close to this, then we might hope to 
reach this solution. 

We will now explicitly give the Ricci-DeTurck flow equations, from which the components 

H 



of R B can be deduced using 



dG ab _ 9R h 
<9A " ~ 2Rab 

dX 

dhjj 

dX 



-2G AC (Rf A - R%Af) 

-2{Rf 5 + R%AfAj - R? A Af - Rf A A A ) (5.5) 



Contracting indices and taking covariant derivatives V with respect to the base metric hij, 
we find, 

=h m Pd m d p G AB + - 

= V*($Gab) -G CD {d l G AD )(d t G C B) + \h k { ™ 1) G CD (d m G CD )(d k G AB ) 
l G BE G AF F Ei 'JF? + i k d k G AB + h^GcDiA™ - A Di )F^d k G AB 



+ h k Z[ ) (A Ci Af ) + A Ci A? - 2A^Af)(d m G CD )(d k G AB ) 
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=h m Pd m d p A<?+ 



-V k Fg + Vi{A c k i k ) + $i($*A° - VP AC) -G AC F B d>G AB + i k Fg 
+\h\ m l) G DE Fgd m G DE - h k ^ r) G DE {A jD - A^)FgPf m 
+\tf? x) {A* D A* + A pD A E - 2A pD A E )Fgd m G DE 
+V« {^h^G DE {Ag - Ag) + G CE (A pD - A pD ) 

+h^ 1) (A mD Ag + A mD A E - 2A mD Ag)(A c - A c ) \ d p G DE ^ 
+Vi(h km l) G D E(A^ E - ti E )(A c k - A c )Fp 



= h m Pdmdphij+.. 

f)h- ' * 1 

= -2R tJ + 2V (i ^) + 2 G AB G CD (d l G CB )(d J G AD ) - G AB F Ak F B 



+ ^\h lk V J (G AB h k ^ 1) d m G AB ) + hikVjih^GABF^A^ - A qA )) 

+ h lk V J (h k L l) (A pA A B + A pA A B - 2A pA A B )d m G AB ) |(m j)\ (5.6) 

where as before the base DeTurck vector field is defined as in (3.4) and we analogously 
define F A = d%A A — djA A and similarly, F t A = d%A A — djA A . In the appendix B to this paper 
we present useful intermediate results that lead to these expressions. 

5.2 Reduced stationary case 

We now make a simple observation, namely that if we require our stationary metric to have 
invariance under the discrete symmetry, 

t^-t, y a ^-y a (5.7) 

then this allows for a consistent truncation of the Harmonic Einstein equation in the sense 
that the Ricci-DeTurck tensor R B U is also invariant. This symmetry implies that all the A A 
vanish, and we explicitly see from the above equations that vanishes, as required for the 
invariance of R B U . Consequently the Ricci-DeTurck flow and Newton methods consistently 
truncate. We term this the 'reduced stationary' case. 

In this case the Ricci-DeTurck tensor considerably simplifies, and has non-zero compo- 
nents, 

Rab = -\v\diGAB) + ^G CD (0 i G AD )(^G C B) - \h km G CD (d m G CD )(d k G AB ) - l ~i k d k G AB 
R% = Rij - V ( 4 - \G AB G CD {d l G CB ){d ] G AD ) - l -h k{i V j) (G AB h km d m G AB ) (5.8) 
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with £ = 0. For a Dirichlet boundary one imposes the induced metric fixing Gab and 
the tangential components of hy. The remaining components of hij are then determined by 
requiring £ J = 0. 

We note that in 4D the 'circularity' theorem [38] implies all stationary vacuum solutions 
may be put in this reduced form, given a condition at a single point which is satisfied for 
asymptotically flat space. More generally all higher dimensional analytic solutions known to 
us are of this form. 

5.3 Boundary conditions for the Killing horizons and axes of symmetry 

We now explicitly give the boundary conditions for the components of our stationary space- 
time metric (5.2) and reference metric (5.4) at the Killing horizons or rotational symmetry 
axes. Recall that for the 4D uniqueness theorems the horizon and symmetry axis play the 
role of boundaries for the Riemannian orbit space, in a very similar manner to our higher di- 
mensional (and cohomogeneity) case here. Therefore we may regard the results in this section 
on the metric behaviour at the horizons and axes as generalising the boundary conditions in 
that context (see for example [6]). They are also consistent with the boundary conditions 
discussed by Harmark using particular coordinates on the base manifold [30] . 

This is an analogous problem to deducing the smoothness condition for a spherically 
symmetric function in spherical polar coordinates at the origin. The function depends only 
on the radial coordinate r, but since it is smooth, meaning in Cartesian coordinates x l it is 
a smooth C°° function of the x l, s, then as r 2 = x l x l it cannot be a smooth function of r, 
but rather is smooth in r 2 . If we require that the function is only C 2 then we deduce that 
the function simply has the Neumann condition df/dr\ r= Q = 0. We may perform a similar 
analysis for a tensor, the only difference being that now the components transform as one 
moves between a chart which manifests smoothness but not the symmetry, and a chart which 
manifests symmetry but not the smoothness. The details of this are straightforward and are 
given in appendix A for a smooth (0, 2) tensor at a Killing horizon or rotational symmetry 
axis. We shall now apply these results to our spacetime metric and reference metric. 

Let us first assume that there is a single Killing horizon, or multiple horizons with common 
normal Killing vector K = T + VL a R a . It is then convenient to change coordinates as, 

t,y a -> t = t, y a = y a -n a t (5.9) 

so that K = d/dt and R a = d/dy a . We note that if R a generates a compact orbit, then the 
coordinate y a is periodic with y a ~ y a + 2tt. Now consider a boundary, either due to the 
vanishing norm of K or a compact R a . We take base coordinates x l = (r,x l ) adapted to the 
boundary so that it lies at r = 0, and decompose the base metric as, 

h ij dx i dx j = Ndr 2 + r N- i drdx i + h Tj dx i dx' j (5.10) 

and likewise for the reference metric where N — > N, iVj —> 2Vj and h~-~- — > h~-~-. 
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Horizon: For a Killing horizon we write the following metric components as, 

G iA = -r 2 f A , A A = rg A , (5.11) 

for A = (t, y a ) and then let X = if A , g A , Gyayb , A A , N , 2Vj , /i^j j be the set of functions 
describing our metric. Let X be the analogous set describing the reference metric. Then the 
results of appendix A imply that for the metric and reference metric to be smooth we require 
the following behaviour; the functions X and X must be smooth functions of r 2 and x l at 
r = 0, and furthermore obey the regularity conditions, 

(f i - K 2 N)\ r=0 = 0, (f i -K 2 N)\ r=0 = (5.12) 

where k is constant and gives the surface gravity with respect to T. We note that since both 
the metric and reference metric are smooth with respect to the same vector field K, which 
is Killing for both, it is the same constant n that must enter the conditions (5.12) above for 
both of them. 

Axis: Consider an axis associated to a vanishing compact R a . Without loss of generality 
choose this to be Rm- Then we choose to write, 

GyN A = r 2 f A , A A = rg A , (5.13) 

and let Y = (f A , g A , G^, G^a , G-^-i , A A , N , N- { , be the set of functions describ- 
ing our metric (where a = 1, . . . , N — 1). Let Y be the set of functions that analogously 
describe the reference metric. Appendix A implies that for a smooth metric and reference 
metric we must have that the metric functions Y and Y are smooth functions of r 2 and x l at 
r = 0, and in addition we require, 

(fy N - N) | r=0 = , {f yN - N) \ r=0 = (5.14) 

Of course we obtain analogous conditions for an axis with respect to a different R a . 

We see that if we have a single Killing vector K = T + Vt a R a normal to all horizons, we 
may use the coordinates (t, y a ) and obtain rather simple boundary conditions. Of course one 
can write these conditions in the original (t, y a ) coordinates. The axis conditions take exactly 
the same form, with replacements t — > t and y a — > y a . If we have multiple Killing horizons 
with different normals, then one must deduce the boundary conditions for each taking coor- 
dinates as in (5.9) with Q a appropriate to each horizon. 

A horizon meeting an axis: It is straightforward to check that the boundary conditions 
at the meeting of a horizon with an axis, or two axes, are compatible with each other. Here 
we will explicitly check this for the metric in the former case, for an axis of Rn, noting that 
one obtains similar results for the other cases and for the reference metric. Take coordinates 
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on the base such that the horizon is at r = and the axis at 9 = 0. Then we write the base 
metric as, 

h ij dx i dx j = Ndr 2 + Md9 2 + r9Adrd9 + r B~ i dvdx i + 9 C- l d9dx i + h Tj dx i dx j (5.15) 

where now i = 1, ...,£) — 3. Then writing, 

Gn = —i" 2 ! i Giyi = r 2 fa , GyNyN = 9 2 g , GyNy& = 6 2 g a G^ N = r 2 9 2 k , 

A A = rp A , A A = rq A , (5.16) 

our arguments from appendix A applied to the horizon r = and to the axis 9 = then 
imply that the set of functions characterising the metric found in the equations (5.15) and 
(5.16) above, N, M, . . . ,k,p A ,q A , together with the remaining components Gyayb and A A , 

must all be smooth functions in r 2 , 9 2 and x % near the meeting point r = 9 = 0. Furthermore 
regularity requires, 

(/ - k 2 N) | r=0 = , and (g-M)\ e=0 = (5.17) 

We see that the conditions from each boundary give rise to a consistent set of behaviours 
above. In particular it implies that two boundaries (a horizon and axis, or two axes) meet in 
the base at right-angles. 

A very important point is that having introduced boundary conditions for the metric we 
must check that the conditions that this implies for the vector £ M are compatible with ensuring 
the elliptic problem in equation (2.6) is well posed with trivial solution. To investigate this 
we must consider our choice of reference metric (5.4), which also is required to be regular and 
hence is subject to the same boundary conditions above for its components. Then one can 
explicitly check that, 

£ r |r=o = 0, d r e\ r=0 = 0, d r £ A | r =o = (5.18) 

both at a horizon and axis of symmetry, which is indeed consistent with well posedness of the 
associated linear problem (2.6) and with a trivial solution. 

One may consider the Ricci-DeTurck flow or Newton method operating on the metric 
g near a horizon. With the choice of reference metric above that has the same isometry K 
and is also regular at the boundaries, the Ricci-DeTurck tensor will be symmetric under K. 
Consider the Ricci-DeTurck tensor in Cartesian coordinates. Since the metric and reference 
metric components will be smooth at the fixed point, then the Ricci-DeTurck tensor will be 
regular there. Thus in our adapted coordinates it will also obey the same regularity conditions 
as the metric and reference metric. In particular, Ricci-DeTurck flow and Newton's method 
will preserve regularity, and therefore we have the very elegant result that they will leave the 
surface gravity constant. The same is true for an axis of symmetry where again regularity is 
preserved. 
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We note that we may view the above smoothness conditions in the weaker sense of 
requiring the Cartesian form to only be C 2 . In this case instead of finding our various 
functions above are smooth in r 2 and x l , we have instead only that these functions obey 
Neumann boundary conditions, ie. vanishing gradient. In addition we also have the regularity 
conditions f\ = k 2 N for horizons, and fyN = N for axes as before. Common folklore is that 
solutions to elliptic problems are analytic, and hence one might expect that only imposing 
such C 2 boundary conditions one would certainly find the stronger smooth behaviour at the 
horizons and axes, and presumably the much stronger analytic behaviour. 

As a final comment we note that for these C°° and C 2 conditions, the elliptic system has 
Neumann boundary conditions for the various functions above, and in addition to this, also 
the constraints, fi = k 2 N for horizons, and fyN = N for axes. One might be concerned that 
these latter conditions should not be imposed in addition to the Neumann conditions, as this 
is 'too much data' for an elliptic problem. However, we emphasize here that this 'fictitious' 
boundary should be viewed as a regular singular point of the equations, due to the singular 
terms arising from the vanishing norm of the Killing vector, and hence a usual counting of 
boundary conditions does not apply. Instead we reiterate that the regularity conditions will 
be preserved by the Ricci-DeTurck flow and Newton method, and thus it is better to think 
of these conditions not as boundary conditions, but rather as a restriction of the problem to 
the class of regular metrics, and using the Ricci-DeTurck flow and Newton method which act 
within that class. 

5.4 An example: Kerr 

It is instructive to consider the simple example of the Kerr solution in light of our discussion 
above. The Kerr solution is in reduced stationary form, so that Af = 0. In the conventional 
Boyer-Lindquist coordinates the Kerr metric takes the form, 

ds 2 = G tt dt 2 + 2G t4> dtd(f) + G H d<p 2 + h rr dr 2 + h ee d9 2 (5.19) 

with fiber metric, 

(A -a 2 sin 2 9) . 2 ((r 2 + a 2 ) 2 - Aa 2 sin 2 9) 
Gtt = ^ , = sin 9 , 

9 (r 2 + a 2 - A) 

G t4> = -a sin 2 6± '- (5.20) 



and base, 



E 
A 



hoe = £ (5.21) 



where the functions A, £ are defined as A = r 2 + a 2 — 2Mr and S = r 2 + a 2 cos 2 9. The 
stationary Killing vector T = and the rotational Killing vector is R = 

The base manifold Ai has coordinates r, 9 upon which the metric components depend 
explicitly. The outer horizon is a boundary of M. and is located at A = where r = = 
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M + yjM 2 - a 2 , and the remaining boundaries are from the axes of rotation at 9 = 0, tt. One 
finds, 

det Gab = ~(a 2 + r(r - 2M)) sin 2 9 (5.22) 

which vanishes at these boundaries. Everywhere in the exterior of the black hole, < r and 
< 9 < tt we have that Gab has Lorentzian signature and hij is Euclidean and smooth. The 
Killing field K = T + QR is tangent to the horizon and timelike near there, where the angular 
velocity of the horizon is given as, Q = q2 " r 2 ■ 

Whilst the 9 coordinate is a regular coordinate on the base manifold at the rotation axes, 
the radial r coordinate is not at the horizon (since A — > there). We therefore define a new 
radial coordinate, p, such that dp = dr/\^A, and p = at the horizon, giving, 

r = M + \/M 2 - a 2 cosh p (5.23) 

so that the components of the base metric are smooth at the horizon boundary. In 
particular in these coordinates the determinant of the base metric, 

h ij dx i dx j = ^dr 2 + ^d9 2 = E {dp 2 + d9 2 ) det = S 2 > r\ (5.24) 

and thus we see that since > the base is indeed a smooth Riemannian manifold everywhere 
on and in the exterior of the horizon. Since we have only one horizon it is convenient to use 
the coordinates discussed above, 

I = t , 4> = 4>-Slt (5.25) 
and then we may confirm that near the horizon, p = 0, we have, 

Gii =-K 2 (h pp \ p=a )p 2 + 0{ P A ) 

G^ = 0(p 2 ), G# = 0{1) 

h pp = h ee = (r\ + a 2 cos 2 9) + 0{p 2 ) 

in accord with our boundary behaviour above, where n is the surface gravity of the Kerr 
solution, 

2 M 2 -a 2 , s 



4M 2 r 2 



At the axis of symmetry 9 = we have, 



G E = 0(1), G u = 0(9 2 ) 
G ii = (h ee \e= )9 2 + O(9 4 ) 



h pp = h ee = (a 2 + (M + V ' M 2 - a 2 cosh/?) 2 ) + 0{9 2 ) 

which again agrees with our calculation of axis boundary behaviour. Likewise the same 
agreement is seen for the axis 9 = tt. 
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6. Example application: 4D rotating black holes in a cavity 



In the above we have set out a numerical framework to find general stationary vacuum black 
holes. We now use the example of AD rotating black holes in a cavity as a toy example to 
demonstrate the methods discussed may be applied straightforwardly in practice. 

We choose to impose Anderson's boundary conditions at the cavity wall [31], where we fix 
the conformal class of the induced metric, and also the trace of the extrinsic curvature of the 
boundary. A canonical choice of cavity is such that the induced metric on the 3 dimensional 
boundary is conformal to the product of time with a round 2-sphere, 

ds\ = -dt 2 + d9 2 + sin 2 9 d(j) 2 (6.1) 

where 4> ~ 4> + 2ir. We must then specify the trace of the extrinsic curvature of the boundary 
which we choose to be constant. If one takes the Schwarzschild solution and cuts it off at 
finite radius, then the extrinsic curvature of the boundary that is introduced is constant 
and positive. Hence in the rotating case we also take this trace of extrinsic curvature to be 
positive. Let us denote this positive constant a. By a global scaling we may choose a to take 
any positive value. For later convenience we choose a = y/2. 

We consider the spacetime to be in the reduced stationary class, with a single Killing 
horizon with spherical topology which is rotating in the <j> direction, so that R = d/d<j) is 
Killing, and K = T + £IR is normal to the horizon. As discussed above it is convenient to use 
the coordinates, 

i=t, (j> = (j)-nt (6.2) 

so that the boundary behaviour at the horizon takes a simple form. Now K = d/dt and 
R = d/d(j). We explicitly write an ansatz for the reduced stationary metric as, 

ds 2 = g^dx^dx" = G A Bdy A dy B + h ij dx i dx j (6.3) 

_ / -r 2 VB r 2 sin 2 9 WB\ a b (a + ^t 2 sin 2 OF 2 r sin 9 f\ . , 
" \r 2 sin 2 9WB sin 2 9SB ) AB { rsin9F B J .. 

where y A = (t, (j>) and x l = (r, 9) and we take the horizon to be located at r = where K is 
null, and the axis of symmetry to be at 9 = 0, tt where R vanishes. Finally we must specify 
the coordinate position of the cavity, and choose this to be at r = 1. Having chosen this the 
base M is then covered by a single chart with domain < r < 1 and < 9 < tt. However, 
we further assume the reflection symmetry 9 — > tt — 9, to reduce the domain to < 9 < tt/2. 

By assumption hij is a Riemannian metric on our coordinate domain. Since K is normal 
to the horizon, we have V > in the neighbourhood of the exterior of the horizon. For 
the metric to be smooth at the horizon, our boundary conditions imply that near r = the 
metric functions {V, S, W, A, B, F} are smooth in r 2 and 9 (or for the metric only being C 2 
we require Neumann boundary conditions in r, and C 2 in 9), with the additional requirement 
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that, 

(VB - k 2 A) | r=0 = (6.4) 

which specifies the surface gravity. At the axis 9 = we require the functions {V, S, W, A, B, F} 
are smooth in 9 2 and r (again for a C 2 metric we require Neumann conditions in 6 and C 2 
in r), with the additional requirement that, 

S\ e=0 = 1 (6.5) 

and likewise for the axis at 9 = ir. The reflection symmetry 9 —> it — 9 imposes that F is 
odd about 9 = ir/2 and the other metric functions are even there. The induced metric at the 
cavity boundary r = 1 in the t, (ft coordinates is conformal to, 

ds 2 B = (-1 + Q 2 sin 2 9)dP + 2SI sin 2 9dtd4> + sin 2 9d4> 2 + d9 2 , (6.6) 

which implies that at r = 1 we have Dirichlet conditions, 

V = (1 -Q 2 sm 2 9), W = n, 5 = 1 (6.7) 

where then the function B gives the conformal factor. The projector onto the induced metric 
on the cavity boundary is = — n^n v where, 

1 f d rsm9F d\ .„ „. 

n =^ [-EZ ( 6 - 8 ) 



sfA \dr B d9 / 

is the outer unit normal to the boundary. The trace of the extrinsic curvature, a, is then 
given as a = L^ v Sj u,n v . The remaining three metric functions B, A, F at r = 1 have boundary 
conditions determined from simultaneously requiring £ r = £ e = together with imposing the 
extrinsic curvature has constant trace a = y/2. These take a somewhat complicated form but 
in essence are coupled oblique boundary conditions. We emphasize that Anderson has proven 
that these give a regular elliptic system [31], ensuring that they give well posed boundary 
conditions for the numerical problem. Indeed we have encountered no problems with the 
cavity boundary in this toy example. 

We see that the physical moduli of a black hole with such boundary conditions, the 
surface gravity k and angular rotation Q, are directly imposed in these boundary conditions. 
We see that for f2 > 1 the horizon Killing field K is no longer timelike. This implies that the 
stationary Killing vector T develops an ergo-region for f2 > 1. We note that for Q < 1 the 
Killing vector T is timelike near the boundary and horizon, although in principle it might 
become spacelike for some intermediate region. However in the results we now present we see 
no evidence of such exotic behaviour. 

For the reference metric we must choose the same form as in (6.3), with the same boundary 
conditions. We make a simple explicit choice, 

J a , UJ v i-r 2 (1 -n 2 r 2 sm 2 9) Sir 2 sin 2 o\ , Aj B (X ()\ , {j , nS 

ds = g uu dx^dx u = v , ' 9 dy A dy B + \ K dx*dx(6.9) 

yM ^ Sir 2 sin 2 9 sin 2 9 J \° V 
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which we note indeed satisfies the smoothness conditions at the horizon and axis, in addition 
to the requirements (6.4) and (6.5) above. The constants k and entering the expression 
above give the surface gravity and angular rotation. One may compute that the trace of the 
extrinsic curvature is equal to k for this reference metric. We note that whilst we will fix the 
trace of the extrinsic curvature for the metric, a = y/2, the same quantity for the reference 
metric is not required to be fixed or equal to a. 

6.1 Numerical results 

We have used Newton's method to solve the Harmonic Einstein equations in this example, 
and have also simulated the Lorentzian stationary Ricci-DeTurck flow. Our aim is not to 
perform high precision numerics, but rather to check that the method behaves as expected. 
We use simple second order finite difference to represent the Harmonic Einstein equations 
and boundary conditions numerically. We have used various resolutions up to 160 x 80 in the 
radial and angular directions respectively, and have checked the convergence of the results is 
consistent with second order scaling. 

Whilst in principle smoothness in r 2 at the horizon should be preserved by Ricci flow 
or Newton's method, numerical accuracy and stability is improved if in addition to taking 
smooth initial data, one also imposes Neumann conditions explicitly on the metric functions 
at the horizon. We take analogous conditions in 9 at the axis boundaries. Likewise whilst the 
conditions (6.4) and (6.5) in principle are again preserved by Ricci flow or Newton's method, 
we also impose these explicitly on the metric functions as boundary conditions to improve 
accuracy. 

We have found solutions using the Newton method starting with the reference metric 
as an initial guess. For all these solutions we confirm that the base metric is smooth and 
Riemannian, although we note that presumably if there was a solution where this was not 
the case we should not expect to have found it. We have checked that the solutions found 
to the Harmonic Einstein equation are indeed Ricci flat, rather than solitons, by ensuring 
that ^ is small and consistent with vanishing in the continuum limit. We note that for a 
reduced stationary metric £ A = 0, and hence the scalar (f> = g^ y ^^, v = hij^ 1 ^ > for hij 
being Riemannian. Hence vanishing <fi implies vanishing Figure 1 shows the maximum 
value of the scalar (j> over our domain for a typical solution, with k = 2 and f2 = 0.6 plotted 
against resolution. We observe similar behaviour for the other values of k and f2. 

Let us define the equatorial radius R eq , horizon area Ah, and cavity boundary area Ab, 



R eq = y/SB\ r=0 e= o , A h = 4tt / d9 \fSBs\nO , A b = 4tt / de^SBsmO 
Let us further define the quantity, 

. = (6.11) 

which we may intuitively think of as giving the ratio of the size of the horizon compared to 
that of the boundary. 



tt/2 

(6. ID) 



=i 
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Figure 1: Plot of the log (base 10) of the maximum value of (f> over our domain for a typical 
solution (Q, = 0.6 and k 2 = 2) against log (base 10) of the number of radial lattice points N. 
The number of angular lattice points is N/2. The figure shows data for resolutions 20 x 10, 
30 x 15, 40 x 20, 60 x 30, 90 x 45, 120 x 60 and 160 x 80. We see that </>, and hence the vector 
field £ are consistent with vanishing in the continuum limit, as we would expect for a Ricci 
flat solution. The same behaviour is seen for all solutions of the Harmonic Einstein equation 
we have obtained. 



We fix Q and then varying k 2 to scan through the moduli space of solutions. We have 
obtained solutions for = 0,0.2,0.4,0.6,0.8,1. and 1.2, although for graphical clarity we 
have not presented all this data in the figures that follow. In figure (2) we plot the quantity 
x against k 2 for fixed £1 = 0.2,0.6,1. Similarly in figure (3) we plot the quantities R eq and 
y 7 A/ Air against k 2 for the same solutions. For Q = the solutions are static, and therefore 
with spherical cavity boundary conditions will simply reproduce the Schwarzschild solution. 
For Schwarzschild one finds, 

n 1 1 x(4-3x) . 



Ax 2 (1 - x) ' eq 2a VT 



x 



and we observe agreement with the data as we expect. We have not plotted the data points 
we have computed for 0, = as these lie on the analytic = curves given in the figures. As 
we turn on the angular rotation, fixing fi, we see deviation from this Schwarzschild behaviour, 
again as we should expect. We note that in the small black hole limit, x — > 0, for fixed f2, we 
expect to recover the Schwarzschild behaviour and we see this is the case. 

For static solutions, so = 0, there is a minimum surface gravity at x = 2/3, and so 
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Figure 2: Plot of x = \J Ah / At against surface gravity k 2 , where Ah is the horizon area. Data 
is shown for numerical solutions with f2 = 0.2,0.6, 1 and 1.2. For f2 = the analytic curve 
is given as a dotted line (the numerical data is omitted as by eye it is indistinguishable from 
this curve). For Q < 1 there is a minimum value of n which divides the branch of solutions 
into small (S) and large (L) black holes. The solutions very close to the minimum are difficult 
to obtain numerically. 

the solutions divide into the small (x < 2/3) and large (x > 2/3) black holes in an analogous 
manner to York's construction where one fixes the induced metric on the cavity wall [39] 
(rather than the conformal class and a). The AdS-Kerr solutions of Carter [40] give an 
indication of what to expect in the rotating case, where in the usual manner we think of AdS 
as being a box. Fixing the AdS length £ = 1, then for £1 < 1 there are small and large solutions 
and a minimum k, and for > 1 there is only one branch of solutions which terminates in 
an extremal solution [37]. In addition AdS-Kerr admits a globally timelike Killing vector for 
Q < 1, and has an ergo-region for Q > 1. 

We recall that for 17 > 1 we must develop an ergo-region for our solutions. Then in 
analogy with AdS-Kerr we might expect that for < £1 < 1 we have small and large black 
holes with a minimum surface gravity, whereas for £1 > 1 there is no minimum surface 
gravity and instead there is an extremal limit where k — > 0. This is indeed borne out by 
our crude numerical results. It would be interesting to confirm with greater accuracy that 
this transition from minimum surface gravity to extremal limit at fixed f2 does indeed occur 
precisely at = 1. Here we have observed that fi = 0.8 has a minimum surface gravity, and 
for fi = 1 we have not found a minimum. We also confirm that for all our solutions with 



- 28 - 




Figure 3: Plots of R eq (red), the equatorial radius, and \J Ah / 47r (blue) against surface gravity 
k 2 . Data is shown for Q = 0.2,0.6 and 1.2, with the analytic curve for Q, = drawn as a 
dotted line. For non zero O and away from the small x limit the horizon becomes deformed 
from the round sphere to become prolate. 

f2 < 1 the vector field K is globally timelike outside the horizon i.e. < 0. 

We briefly comment on finding the solutions using Newton's method. We begin by finding 
a solution for some k, starting with the reference metric as an initial guess. For SI < 0.6 
starting with k 2 = 2 this yields a small black hole solution with boundary trace of extrinsic 
curvature a = v2 as required. We then move along the branch of solutions by perturbing 
the solution and reference metric to yield a good approximation to a solution with nearby n 
and the same O and a. Using this method we may quickly scan along a branch of solutions 
to the minimum value of k. We could try to find an initial guess in the basin of attraction 
of a large solution and then scan along this branch. However the approach we take is to 
extrapolate our small black hole solutions near the minimum k to gain a good guess for a 
large solution, with k just greater than the minimum. Using this one finds a large solution, 
and can then scan along the large branch. The Newton method struggles to find solutions 
very close to the minimum value of n. At the minimum there is a normalisable zero mode 
of the linearised Harmonic Einstein equation, and near to it there is a low lying mode that 
renders the linear operator that must be inverted in Newton's method rather ill conditioned. 
This is why in figures 2 and 3 for SI = 0.2 and 0.6 the section of the curves connecting the 
small and large branches are missing. In order to find black holes with Q > 0.6 we have found 
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that extrapolating these solutions together with their reference metrics for some fixed k, say 
k 2 = 2.5, on the small black hole branch then allows an initial guess for a higher value of O 
to be found. Having found this new solution branch, one can proceed to iterate this method 
to find larger f2 solutions. 




(a) O = 0: Small branch (b) O = 0: Large branch 

Figure 4: Figure showing the upper quadrant of a polar section of the embedding of the 
horizon (red) and cavity boundary (blue) as a function of k 2 for fixed $7 = 0. A slice at 
constant n should be rotated about the vertical axis, and then reflected in the horizontal 
plane to obtain the surface of revolution of the embedding in M 3 . The solutions plotted 
are for the small black holes (left frame) and the large black holes (right) frame, and being 
Schwarzschild (as = 0) the horizon and cavity boundary are spherical. This figure should 
be contrasted with the later figures which give the embeddings for Vt > 0. 

We may isometrically embed the horizon and cavity boundary of these solutions into the 
Euclidean space M 3 as surfaces of revolution. Taking a polar slice through these embeddings 
then gives curves in two dimensions representing the geometry of the horizon and boundary. 
In figures 4, 5 and 6 we plot the upper quadrant (i.e. < 9 < tt/2) of these two dimensional 
embedding curves against k 2 for Q = 0,0.6 and 1. Since for O = the solutions are simply 
Schwarzschild, figure 4 is included simply for comparison. For the solutions found with < 1 
there is both a small and a large branch. Interestingly for O = 0.6 we see that for the large 
solutions with sufficient k, the geometry of the horizon is such that its vertical extent in the 
embedding is actually greater than that of the cavity indicating the geometry of the solution 
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2.0 



2.0 



(a) fl = 0.6: Small branch (b) £7 = 0.6: Large branch 

Figure 5: This figure shows the horizon and cavity boundary embeddings of the small rotating 
black holes (left frame) and large rotating black holes (right frame) with f2 = 0.6, and should 
be contrasted with the previous figure depicting Schwarzschild. Interestingly we find for the 
large black holes the vertical size of the embedding of the horizon may exceed the vertical 
size of the cavity embedding. 

is rather exotic. For = 1 we find that the horizon may not be embedded isometrically into 
M 3 for k? < 0.7 in an analogous manner to the Kerr solution sufficiently near extremality [41]. 

We conclude our discussion of these solutions by studying certain stationary Ricci- 
DeTurck flows. From previous work [22, 17] we know that for static black holes in a spherical 
cavity with fixed induced metric the small black hole is always unstable having a Euclidean 
negative mode [35, 42] whereas the large black hole is stable. There are two flows for a small 
black hole generated by this unstable mode. Perturbing the small black hole by the negative 
mode in one sense generates a flow that asymptotes to the large black hole. By reversing 
the sign of this perturbation one then generates a second flow that at finite time develops a 
singularity where the horizon shrinks to zero size. In [22] it was argued that by an appropriate 
surgery on the manifold one may continue the flow to flat space. 

We note that whilst for fixed induced metric the static spherically symmetric problem 
(continued to Euclidean signature) is well posed [23], in our cohomogeneity two example we 
must take Anderson's boundary conditions to obtain a regular elliptic system. Interestingly 
we find that for these boundary conditions and fi < 1 the large black holes now have one 
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Figure 6: Figure showing the horizon and cavity boundary embeddings for rotating black 
holes with $7 = 1. This branch actually extends to lower k although the horizon cannot be 
globally embedded in M 3 isometrically for k 2 < 0.7. 



mode of instability for Ricci-DeTurck flow, and the small black holes tested have two. We 
emphasize that since we are considering the flow in the space of stationary Lorentzian metrics 
these are not Euclidean negative modes, but are rather negative eigenvalue eigenmodes of the 
Lichnerowicz operator restricted to our class of stationary spacetimes. This implies that 
to find the large black holes one must tune a one parameter family of initial geometries to 
reach these unstable fixed points. For the small black holes presumably one must tune two 
parameters although we have not tried this. 

We now give an example of the Ricci flow algorithm to find a large black hole by tuning 
a one parameter family of data. We present results for Q = 0.6 and k 2 = 2, although we 
expect to see qualitatively similar results for other Q < 1 and n greater than the minimum 
surface gravity. Let us take a one parameter family of initial data to be, 



-r 



[1 -n 2 r 2 sin 2 9) (1 + 5 f) Clr 2 sin 2 e\ A B (\{l + 8f) 



ds = v '9 9,/ 9. dy A dy B +\^ K J ' ",2 dx'da? 

\ il,r 2 sin 2 9 sin 2 6 yy I 1 + 6 (rf) 2 ' 



AB \ / ij 



with / = (l — r 2 ) and where 5 is the parameter we may tune. We choose k 2 = 2 and Q = 0.6, 
and then the trace of the extrinsic curvature a = y/2 as required, independent of the value 
of 5. Likewise the constants n and Q entering the expression above give the surface gravity 
and angular velocity independent of 5. The above metrics also satisfy all our other remaining 
boundary requirements discussed above. The reference metric is as before in equation (6.9). 
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(c) S = 2.9 



Figure 7: Figure showing the embeddings of the horizon (red) and cavity boundary (blue) 
against Ricci flow time A for three values of the parameter S. For a critical value 5 ~ 2.9 the 
flow tends to the large rotating black hole solution with ft = 0.6 and k 2 = 2. For smaller 8 
the horizon and boundary appear to collapse to a singularity at finite flow time. Conversely 
for greater 5 they appear to expand in an unbounded manner under Ricci flow. 

We find that we may reach the critical point for 5 ~ 2.9. In figure 7 we depict embeddings 
of the horizon and boundary for 3 values of 5, one for slightly less than the critical value 
(5 = 2.8), one for slightly greater (5 = 3.0), and one close to the critical value (5 = 2.9). 
We see that for 5 smaller than the critical value the horizon and cavity boundary appear to 
collapse to a singularity at finite flow time. For 5 greater than the critical value the horizon 
and boundary appear to expand forever under the Ricci flow. For an appropriately tuned 
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flow we may approach the large rotating black hole solution with Q, = 0.6 and k 2 = 2. This 
is the same solution found by the Newton method. Indeed we see the horizon and cavity 
embeddings in the tuned flow in this figure match those presented in figure 5 for k? = 2. 

7. Summary 

In this paper we have discussed applying the Harmonic Einstein equation with DeTurck term 
to the case of Lorentzian static and stationary spacetimes, where Ricci flow or Newton's 
method acting on these spacetimes can be used as algorithms to solve this equation and yield 
solutions to the static or stationary vacuum Einstein equations. 

Previously static spacetimes have been treated by continuing to periodic Euclidean time, 
which can remove boundaries associated to non-extremal horizons [17]. The Harmonic Ein- 
stein equation is then solved as an elliptic boundary value problem on a Riemannian manifold, 
for example, using Ricci flow or Newton's method. 

We have firstly shown how to view this procedure for static spacetimes from an entirely 
Lorentzian point of view. It is straightforward to see, with a suitable choice of reference metric, 
that the Lorentzian Harmonic Einstein equation consistently truncates to the static case and 
gives an elliptic system. Whilst the horizon must be treated as a boundary now in Lorentzian 
signature, by considering the Euclidean continuation, it is easy to deduce Lorentzian boundary 
conditions for regularity at the horizon. Ricci flow on the space of static spacetimes, or 
Newton's method can then be applied to solve for static solutions of the vacuum Einstein 
equations. 

Having treated the static case from a manifestly Lorentzian perspective we then con- 
sidered the stationary case for globally timelike stationary Killing vector. As in the static 
Lorentzian case it is straightforward to see that for a suitable choice of reference metric, the 
Harmonic Einstein equation consistently truncates to such stationary spacetimes, again giv- 
ing an elliptic system. However, we learn that in general, if the stationary Killing field fails 
to be timelike, such as at horizons or in ergo-regions, it is unclear the system will remain 
elliptic. 

Intuitively we may say that the challenge to ellipticity in the Lorentzian setting is spatial 
gradients in the direction of motion of the horizon. However for a spacetime to be stationary 
there can be no such gradients as is formalised in the Rigidity theorems. If there were, radi- 
ation would be emitted in contradiction to stationarity. Motivated by Rigidity we consider a 
broad class of stationary spacetimes which are given by a fibration of the orbits of the sta- 
tionary Killing vector, together with orbits of other mutually commuting Killing vectors that 
generate rotation or translation and also commute with the stationary field. We then require 
that the Killing vector normal to a rotating horizon is some linear combination of these vector 
fields. This is the same class considered by Harmark in order to classify higher dimensional 
black holes [30]. Further motivated by the stationary uniqueness theorems we make the key 
assumption that the base of this fibration is a Riemannian manifold with boundaries at the 
horizon and any axes of the commuting rotational symmetries. It is then straightforward to 



- 34 - 



see that for a reference metric of the same form, the Harmonic Einstein equation consistently 
truncates to this class of stationary spacetimes. We provide the necessary boundary condi- 
tions at the Killing horizons or axes of symmetry and demonstrate they are consistent with 
obtaining vacuum solutions of the Einstein equations rather than only Ricci solitons. The 
Ricci-DeTurck flow is then parabolic on this class of Lorentzian stationary spacetimes, and 
gives an explicit algorithm to solve it, as does Newton's method. 

In order to demonstrate that these methods may be applied in practice we have considered 
a very simple example where analytic solutions are not known, namely a 4D rotating black 
hole in a spherical cavity. We have used Anderson's boundary conditions at the cavity wall, 
noting that the Dirichlet conditions (i.e. fixed induced metric) often discussed in the static 
spherically symmetric case are not well posed. Solutions were found using Newton's method. 
Ricci flows of these Lorentzian stationary black hole spacetimes were also performed, and their 
use demonstrated to construct solutions. We believe these are the first example of Lorentzian 
Ricci flows that have been studied. 

Of course this toy example is only of cohomogeneity two and is 4 dimensional. The 
main purpose of this paper is precisely to give an elegant geometric framework that applies 
beyond these cases, and it will be interesting in future work to test these methods in the 
many physically interesting situations where stationary black holes have been conjectured to 
exist. 
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A. Boundary conditions for horizons and axes 

Our elliptic problem is posed on a base M. which has boundaries associated to Killing horizons 
and axes of symmetry. We must explicitly give boundary conditions that ensure regularity 
for these. In order to derive the necessary conditions, we perform a similar calculation to 
that in section 3 where we used the smooth Euclidean continuation of a static black hole to 
deduce boundary conditions for the Lorentzian problem. 

We are interested primarily in regularity of the metric. However we have argued that 
in the static Lorentzian case, as is clear from the Euclidean picture, that the Ricci-DeTurck 
tensor also shares the same regularity. Therefore we perform the analysis below for a general 
(0, 2) tensor J which one can take to be the metric or Ricci-DeTurck tensor. 
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The procedure is best illustrated by a simple example. Consider a smooth (0, 2) tensor, 
J, which is symmetric with respect to a vector field R that generates U(l) orbits with period 
2-7T and has fixed action at some point p. We will then derive regularity conditions for the 
components of J in a chart with 'polar' coordinates (r, a) adapted to the symmetry so that 
R = g^, and hence the components do not depend on a. 

We take a to have period 2tt, and the fixed point of R is r = 0. To proceed the tensor J 
is first written in 'Cartesian coordinates' (a, b) that do not manifest the U(l) symmetry, but 
where J has components which are C°° smooth everywhere, including the fixed point; 

J = N(a, b)da 2 + M(a, b)db 2 + K(a, b)dadb (A.l) 
We then introduce the 'polar coordinates' r,a that make explicit the U(l) symmetry as, 

a = rsina, 6 = rcosa (A. 2) 

Let us write J in this polar chart as, 

ds 2 = r 2 A(r)da 2 + B{r)dr 2 + r 3 C(r)dr da (A. 3) 

where the metric functions A, B, C now depend only on v. The symmetry conditions, — 
§^ = = then translate into conditions on the original metric functions N, M and K 
and evaluation of these gives the regularity conditions, N = M and K = at the fixed point 
a = b = 0. When expressed in terms of the polar metric functions these conditions become, 
[A — B)\ r= Q = 0. Consider the function A which can be given in terms of the Cartesian 
components as, 

A = a2 l +h2 (a 2 M + b 2 N - abK) (A.4) 

As a consequence of the above regularity conditions, we see that A is a smooth function of 
a, b. Since it is only a function of r, and r 2 = a 2 + b 2 , it follows that A is a smooth function 
of r 2 . Similarly one finds B, C are also smooth in r 2 . 

Thus we conclude that if we write the tensor J is the polar chart as in (A. 3), we may 
think of r = as a boundary, where we impose that A,B,C are smooth in r 2 there, and 
(A-B)\ r=0 = 0. 

A.l Regularity and smoothness at a Killing horizon 

A Killing horizon implies the existence of a normal Killing field K, whose isometry group is 
K, with a fixed point at the bifurcation surface, and whose orbits close on the future and past 
horizons. We consider a smooth (0, 2) tensor that is symmetric under K, in a chart adapted 
to the symmetry which covers the exterior of the Killing horizon. The fixed point may be 
regarded as a boundary of this chart, and we determine regularity conditions on the tensor 
components there. 
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We begin in smooth Cartesian coordinates with a (0, 2) tensor, J, written in components 

as, 

J = Nda 2 + Mdb 2 + Udadb + Q i dadx i + R i dbdx i + T ij dx i dx j (A.5) 

where i = 1, — 2. We take a Killing horizon with respect to the Killing vector i*T to be 
located at a = b and a = —b with bifurcation surface a = b = 0. Since J is smooth at the 
horizon these component functions are C°° in the neighbourhood of the horizon. The horizon 
Killing symmetry is not manifest in these coordinates and in analogy with the toy example, 
we now change to hyperbolic coordinates, 

a = rsinliKi, b = rcosh Kt (A. 6) 

so that K = d/dt and r = is the bifurcation surface, with k a constant related to the 
normalization of K and giving the surface gravity. We write the metric in this polar form as, 

J = -r 2 Adt 2 + Bdr 2 + r 3 Cdrdt + rFidrdx 1 + r 2 G i dtdx i + T ij dx i dx j (A.7) 

where the component functions are independent of t. Repeating the analysis outlined in 
the toy example one arrives at the conclusion that the functions A, B,C, Fi,Gi,Tij depend 
smoothly on r 2 and x l , together with the regularity condition, 

A\ r=0 = K 2 B\ r=0 (A.8) 

Thus we see explicitly that the regularity in the chart A.7 depends on the normalization of 
K, and hence the surface gravity. If we take the tensor J to be the metric, we deduce the 
regularity conditions on the metric at the Killing horizon. Taking J to be the Ricci-DeTurck 
tensor we see the behaviour it will exhibit if it shares the symmetry and is regular. 

A. 2 Axis of rotation 

We now consider a (0, 2) tensor which is symmetric under a Killing field R which generates 
rotation about an axis with period 2ir. In a chart which manifests the symmetry the axis is 
fixed under the U(l) action, and may be regarded as a boundary for the chart. We determine 
the regularity conditions for the components in this chart there. This case is very close to 
the toy example before. 

We begin with a Cartesian line element of the form 

J = Nda 2 + Mdb 2 + Udadb + Q i dadx i + R i dbdx i + T ij dx i dx j (A. 9) 

and the component functions depend smoothly on a, b, x l in the neighbourhood of the axis 
which we take to be a = b = 0. We now change to polar coordinates defined by, 



a = r sin a , b = r cos a 



(A.10) 



where a has period 2ir and R = d/da and r = is the axis. In these coordinates we write 
the tensor as, 

J = r 2 Ada 2 + Bdr 2 + r 3 Cdrda + rF i drdx i + r 2 G i dadx i + T ij dx i dx j (A.ll) 

and the symmetry is manifest so the metric functions are independent of a. Repeating the 
analysis outlined in the toy example, one finds the metric functions A,B,C,Fi,Gi,Tij are 
smooth functions of r 2 and x l , together with the regularity condition, 

A\ r=0 = B\ r=0 (A.12) 
B. Connection Components and Flow Equations 

In this appendix, we give the connection components of the metric 5.2 together with the 
components of the Ricci tensor and £ vector. The Christoffel symbols are given by, 

T) k = t) k + \h im A Aj F A m + h im A Ak Ff m - \w™A A Ap m G AC 

Tab = 

ri = - l -A A iG BC F% + ^AfdjGBc + lG AC diG BC 

T A = -A A f% + l -A Ak A Bl F^ 3 + l -A Ak A B] F B + d u A A + G AB A%d 3) G BC + l -A Ak A B A? d k G DB 
r ec* = -^A At diG B c 

r}A = -\h lk Af d k G AB + l -h ik G AB F B (B.l) 

where Y l - k is the Christoffel connection of the 'sub metric' hij and F-j = d{A A — djAf = 
ViA A — V jA A . The covariant derivative in the latter equation, Vj, is defined with respect to 
the connection T l - k of hij (and is metric compatible with respect to hij). Using these results 
one finds for the decomposition of the Ricci tensor, 

Rab = \wViidjGAB) - ^G CD h ip (d p G AB )(diGcD) 

+ \h ij G CD {djG CB )(diG AD ) + ^h mi h^G BE G AF F^F[ p 

RiA - R AB Af = ^h jk G AB VjFg + h? k Fp k G AB + ^h jTn G CD G AB F? m djGcD 

Rij + R AB A A Af - R Ai A A - R A3 A A = R i3 - ^G CB V 3 (diG C B) 

+ hG CD G BA {d ] G DA ){d l GcB) + \h km G AB F A m F B 
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where Rij is the Ricci tensor computed with respect to f* fc . The DeTurck vector £ M = 
g Xu (T^ u — T^ u ) decomposes as, 

e = i k - \G AB h km d m G AB + \G AB h k ^ x) d m G AB + Wh km x) G AB {A A - Af)Ff m 

+±h i ih k ™ ) (AfAf + AfAf - 2AfAf)d m G AB 
f = l G AB A Cj d . GAB _ \q ab A c i djG AB + WtfiAf - <7 t Af) 

A h ^A Ck (AfAf + AfAf - 2A A Af)d k G AB + h^A Ck G AB (A A - A A )F B 
+h^G CB {A A — Af)djG AB (B.2) 

where £ fc = h l] (T k j — T k -) and as usual, an overbar indicates that the quantity in question is 
evaluated in the reference metric. 

As discussed in the main text, the flow equations for the various metric components of 
interest decompose as, 

-2R AB + 2V(^£b) 

-2G AC (R jA - R AB Af) + 2G AC (V (A ^ - V (A £ B) A B ) 

-2{Rij + R AB A A Af - R lA A A - R ]A A A ) 
+2(V { ^) + V (A £ B) A A A B - V {B ^Af - V^Af) 

This form is particularly useful as the linear combinations of the components of V^,,) 
that arise take a relatively simple form. Explicitly one finds that, 

2V (a £b) = i k d k G AB - l 2 G CD (d k G CD ){d k G AB ) + 2 G CD h k ^ 1) (d m G C D)(d k G AB ) 
+ h^GcDF^ - A> c )d k G AB 
+ ^h^A^Af + A* C A? - 2A lC Af)(d m Gc D )(d k G AB ) 



dG AB 
dX 

dX 
dhjj 

l)x 
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2(V (i d) - V {A iB)Af) = G AC V t (A^i k ) + G AC %(V P A^ - V P AC) 
+ GAcVi ^-h^ 1)G DE {A c m - Ag) + G CE {A? D - A pD ) 

+±h kp 1) (A™ D Ag + A mD A E l - 2A mD Ag)(A% - A c k )^ d p G DE ^ 
+ GAcV^h^GoEiA^ - A^ E ){A c k - A^Fp 
+ G A d k Fg - l 2 G AC G DE Fg{d k G DE ) 

+ 2 G acG de F ki d m G de 

- hp 1} G AC G DE (A jD - ti D )FgF E m 

+ h k ^ l) G AC (AP D A E + A pD A E - 2A pD A E )Fgd m G DE 

2(V (i Ci) + V {A £ B) A?Af - V [B ^Af - V (B ^Af) = 2V ( 4 + G AB G CD (d l G CB ){d J G AD ) 

- G AB V t (d 3 G AB ) + (^hikV^G^h^drnGAB) + h ik Vj (h^G AB F B m {A^ A - A* A )) 

+\h ik V j {h\™ x) {Av A A B + A* A A* - 2A? A A B )d m G AB ) + (it* j)) 

where we note that in these latter three expressions, all 'A term' base indices have been 
contracted with the base metric h, L j as appropriate. Using these results, one arrives at the 
flow equations in the main body of the paper, contracted in the same manner. 
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